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We apply the relativistic chiral Lagrangian to the nuclear equation of state. An ef- 
fective chiral power expansion scheme, which is constructed to work around nuclear 
saturation density, is presented. The leading and subleading terms are evaluated 
and are shown to provide an equation of state with excellent saturation properties. 
Our saturation mechanism is found to probe details of the nuclear pion dynamics. 



1 Introduction 



^ I The nuclear equation of state awaits a fundamental theoretical description. 

in ■ It is therefore important to develop the appropriate form of effective chiral 

perturbation theory (ExPT), a promising tool of modern nuclear physics, and 
vQ . apply it to the nuclear many body problem systematically. 

f^ I The key element of any microscopic theory for the nuclear equation of 

On ■ state is the elementary nucleon-nucleon scattering process. In the context of 

On I chiral perturbation theory this problem was first addressed by Weinberg who 

proposed to dfjdve a-ciiiral nucleon-nucleon potential in time ordered pertur- 
bation theoryld'cl. Inl3'cl we proposed an alternative scheme with chiral power 
counting rules applied directly to the nucleon-nucleon scattering amplitude. 
^ ' This approach applies the manifest covariant version of the chiral Lagrangian 

P3 . and relies on the crucial observation that the chiral power counting rules can 

be generalized for 2-nucleon reducible diagrams. Non-perturbative effects like 
the pseudo bound state pole in the ^Sq channel are generated by properly 
renormalized local two-nucleon vertices. A similar approach was suggested by 
Kaplan, Savage and WiseO in the context of the chiral potential approach. 
Unfortunately the straightforward application of the extended counting rules 
to the s-wave nucleon-nucleon scattering phase shifts is hampered by a poor 
convergence and inconsistency with empirical dataj. In particular, unitarity is 
violated strongly even at small energy where one would expect this scheme to 
work welQ. This is, at first glance, a surprising result since it is long known 
that for the one pion exchange Yukawa potential pertuihation theory works 
fine if the coupling strength is put to its empirical valueu'ErLi. In fact in the ^Sq 
channel, the failure of a naive expansion is linked to the presence of a second 
strong subthreshold singularity reflecting a subtle interplay of intermediate 
attraction and short range repulsiorcl. A properly modified unitary implemen- 
tation of the extended chiral counting rules describes empirical s-wave elastic 



and inelastic nucleon-p.ucleon scattering in the spin singlet channel accurately 
up to E^iab ^ 600 Me\H. 

We point out that counting rules for the vacuum nucleon-nucleon scatter- 
ing amplitude greatly facilitate ±he construction of chiral power counting rules 
for the nuclear matter problenxHI. Basically the density expansion or equiva- 
lently the multiple scattering expansion can be combined smoothly with the 
chiral expansion where one identifies fc^? ~ Q as a further small scale. Such 
an identification leads to a systematic resummation of the density expansion 
avoiding the expansion in large ratios like kp/mT^ but exploiting the chiral 
mass gap and expand in small ratios like kp/mx or m-„/my^ with rrix — ttin- 

2 Primer on relativistic chiral nucleon-nucleon scattering 

In this section we recall basic elements required for a successful application 
of the relativistic chiral Lagrangian to the baryon sector. The systematic 
construction of the Lagrangian density in accordance with chiral symmetcy 
was presented in great detail by Weinberg using stereographic coordinatesll^. 
There are two problems when evaluating the relativistic chiral Lagrangian in 
the baryon sector. First, any covariant derivative acting on the nucleon field 
produces the large nucleon mass and therefore must be assigned the minimal 
chiral power zero. Thus an infinite. tmr.er of interaction terms needs to be 
evaluated at given finite chiral orderOJH. Second, the straightforward evalu- 
ation of relativistic diagrams involving nucleon propagators generates positive 
powers of the large nucleon mass from loop momenta larger than the nucleon 
mass. The chiral power counting rules are spoiled. The relativistic scheme 
appeared impractical and the heavy mass formulation of xPT, which over- 
comes both problems by performing a non relativistic 1/m expansion at the 
level of the Lagrangiap_jdensity, was developed Eil and applied extensively in 
the one- nucleon sector E3. 

We point out that in the two-nucleon channel it will 
work with the relativistic form of the chiral Lagrangian 
that in the heavy baryon formulation two-particle reducible diagrams are ill 
defined. A second reason is that a proper treatment of the pion production 
process pp — *■ ppn^ requires relativistic kinematica. Note that an alterna- 
tive scheme to tame potentially ill defined two-particle reducible diagrams is 
the use of time-ordered perturbation theory as suggested by WeinberjI. Since 
covariance is not manifest in time-ordered perturbation theory we prefer the 
relativistic scheme. Here we outline how the two aforementioned problems can 
be overcome within the framework of the relativistic chiral Lagrangian. Of 
course we expect the relativistic scheme to be equivalent to the more familiar 



B advantageous to 
. One reason is 



heavy fcrmion formulation of chiral perturbation theory in the one-nucleon sec- 
tor. Rather than performing the l/m expansion at the level of the Lagrangian 
density we suggest to work out this expansion explicitly at the level of indi- 
vidual relativistic Feynnian diagrams. Therewith we avoid the heavy baryon 
chiral Lagrangian with its known artifacts in the two-nucleon sector. Here 
it is convenient to consider the chiral Lagrangian (at an intermediate stage) 
to represent a finite cutoff theory. The finite cutoff, A <C m, is required to 
restrict the nucleon virtuality inside a given loop diagram such as to justify 
the \/m expansion. Obviously a finite cutoff has to be introduced with great 
care not to break any chiral Ward identity. Note that one may alternatively 
apply dimensional rcgularization. We emphasize that obviously if dimensional 
regularization is applied one has to first expand in l/m and then perform the 
loop integration. The two steps do not commute here. 

We need to outline how to arrange higher order terms of the relativistic 
chiral Lagrangian in the nucleon sector. To have a practical scheme one would 
like to have only a finite number of interaction terms with a given number of 
'small' derivatives involved. The known problem of the 'large' nucleon mass 
generated by time derivatives acting on a nucleon field has to be taken care of. 
This problem is cured easily as follows. Consider a covariant derivative ^^{x) 
acting on a nucleon field N{x). We need to construct a 'counter interaction' 
which leads to a small renormalized vertex. This rearrangement is illustrated 
at hand of the simple case when the Lorcntz index fi is saturated by another 
covariant derivative acting on the same nucleon field. In this case the counter 
interaction is readily found 

V^'{x)Vf,{x)N{x) ^ {V^'{x)V^{x)+m'^) N{x) 

^ - (p2 _ m^) ^ Q (1) 

and the 'renormalized' vertex is proportional to the virtuality of the nucleon 
4- momentum which is of 'minimal' order Q. It is important to observe that 
this counter interaction must have occurred already in the chiral Lagrangian 
to some lower order since it involves fewer derivatives. Note that a similar 
construction works if two covariant derivatives act on different nucleon fields. 
We are left to consider the case in which the Lorcntz index /i of the co- 
variant derivative I?^ is contracted either with the Lorcntz index of a Dirac 
matrix or the Lorentz index of a derivative acting on a pion field. Interaction 
vertices in which a covariant nucleon derivative V^ is saturated with a pionic 
four vector D'^ or with the Lorentz index of its covariant derivative behave well 
since they are already suppressed by the small derivative acting on the pion 
field. Here one has to assign the nucleon derivative the chiral power Q^. The 
large nucleon mass present is not hazardous provided that one may identify the 



natural scale A of chiral interactions with the nuclcon mass. With A ~ mp this 
seems indeed reasonable. The remaining case in which the nucleon derivative 
couples to the Lorentz index of a Dirac matrix does not generate an infinite 
tower of interaction terms since at a given number of nucleon fields, n^, the 
number of available Lorentz indices provided by Dirac matrices is at most Ui . 

The above arguments show how to organize the infinite tower of chiral 
interaction terms such that there is always only a finite number of terms with a 
given number of 'small' derivatives. Thus at any given chiral order only a finite 
number of diagrams need to be evaluated. The result of this reorganization 
of the relativistic chiral Lagrangian leads to the applicability of Weinberg's 
counting rules provided that individual diagrams are expanded properly. 

In the two-nucleon sector the chiral Lagrangian leads to a well defined 
expansion scheme for the Bethe-Salpeter kernel K. For the unitary iteration of 
K, as induced by the Bethe-Salpeter equation, Weinberg's counting rules are 
not applicable due to the s-channel unitarity cut. Nevertheless a generalized 
counting for two-particle reducible diagrams can be established. For a given 
diagram each pair of intermediate nucleons causes a reduction of one chiral 
power as compared to the 'naive' chiral poweru. The argument goes as follows: 
consider the once iterated Bethe-Salpeter kernel 

KGK^f^^ '-^^ '-^^^ (2) 

with W^ = s. The object G represents the two-nucleon propagator and K the 
Bethe-Salpeter kernel evaluated according to Weinberg's counting rules to a 
given chiral order. It is instructive first to study the s-channcl spectral density, 
p{s), of this contribution. Since by assumption the Bethe-Salpeter kernel K 
is 2-nucleon irreducible the s-channel spectral density picks up strength only 
from the pinch singularity generated by the two nucleon propagators in (0). 
In the center of mass frame with W — Q one finds: 

p{s) = <iKGK^ I f- ffiK{l)K{l) (3) 

with|^ s — 4: (r?!^ + P^) • Expression ^ can now be used to introduce the chiral 

"Of course the iterated Bethe-Salpeter kernel does not necessarily satisfy a s-channel dis- 
persion relation. However, we point out that only the part which provides strength for the 
s-channel spectral density requires special attention and modified power counting rules. Ob- 
viously a part which does not generate strength for the s-channel spectral density has no 
pinch singularity leading to a s-channel cut so that standard counting rules apply. 



power, I/, of the s-channel spectral density, p(s), 

v=\ + 2vK (4) 

with the chiral power vk of the Bethc-Salpeter kernel K as given by Weinberg's 
power counting rules. In (H) we exploit the important observation that the 
small scale momentum p ^ Q occurs in the spectral density. 

To have a useful scheme we also need to derive counting rules for the real 
part of the scattering amplitude. Here we evoke causality which relates the 
real part of a given Feynman diagram to its imaginary part by means of a 
dispersion relation 

I{s)^^KGK = -V r ds—^ . (5) 

T^ JAm'^ S-S-ie 

Of course we must not assume the dispersion relation dq) to be finite and well 
defined. It may be ultraviolet divergentn. However, we recall that we are free 
to consider our chiral Lagrangian as a finite cutoff theory. So let us introduce 
a covariant cutoff A into the dispersion relation (P). If the cutoff A is to be 
interpreted as a bound for the maximal virtuality of nucleons one should first 
remove the nucleon rest mass. Therefore we take the cutoff A to restrict the 
momentum p < A rather than the Mandelstam variable s: 

A 2 

i{p,A)^-rf dp^4^. (6) 

The natural guess for the chiral power of the integral /(p) would be the chiral 
power of its imaginary part. However, the integral (o) as it stands does not 
permit this chiral power yet. Though the spectral density p(p) exhibits only 
small scales as shown above, the cutoff scale A prevents the conclusion that 
the integral behaves like a small scale to the power \^2vk- We point out 
that the counting rule (13) can nevertheless be applied to the real part of 
the loop function, however, only after an appropriate number of subtractions 
at p^ = — /i^ < with ^1 ^ ttIt^. The subtraction is required to render the 
dispersion integral finite in the large cutoff limit: 

TT Jo P'^ -P'^ VP +M / 



''Note that for a given contribution one has to specify whether the dispersion relation holds 
at fixed Mandelstam variable t or u. This technical detail, however, does not affect our 
argument. 



The subtraction polynomial can always be absorbed into the local two-nucleon 
vertices. Note that the performed subtractions do not necessarily simplify our 
attempt to apply dimensional counting rules since this procedure introduces 
a new scale, /x, the subtraction point. Nevertheless we now can exploit the 
freedom to choose the subtraction point. If we insist on a 'small' subtraction 
scale n ~ m^r of the order of the pion mass one can apply dimensional power 
counting since the integral is finite in the large cutoff limit by construction. 
Thus we arrive at our desired chiral counting rule (13) for 2-particle reducible 
diagrams: each pair of nucleon propagators which exhibits a s-channel unitarity 
cut gets a chiral power —3 in contrast to the 'naive' power —2. We stress, 
first, that this counting rule holds only if one introduces a 'small' subtraction 
scale fi '^ niTr at the level of the Bethe-Salpeter equation and second, that 
the subtraction scale is independent of the intrinsic cutoff. We will refer to 
this counting rule as the 'L'-counting rule since for the one pion exchange 
interaction the chiral power of the n-th iterated interaction is simply given by 
the number of loops. InQ the reader may find explicit examples which confirm 
our L-counting rule (0). 

The non-perturbative structures like the pseudo bound state in the spin 
singlet channel are generated naturally since the bare 2-nucleon vertex is renor- 
malized strongly so that it effectively carries chiral power minus one. This can 
be seen as follows. Consider the leading order interaction 

J^iN{x) = i 5 {N{x) 75 C TT2 N\x)) {N\x) t^ t C^Ss N{x)) . (8) 

with the charge conjugation matrix C — i 70 72 and the isospin Pauli matri- 
ces Ti . At tree level the coupling g can be interpreted in terms of the s-wave 
scattering lengths a(^5o) with 2mg = —Air aQ^S^). We assume here that the 
energy dependent 4-nucleon bare coupling strengths, g(jp), is natural. Accord- 
ing to our subtraction scheme, which is required for the manifestation of the 
L-counting rule, the bare coupling g is renormalized by the one loop bubble J 
at leading order: 

g-^p") - J{p\K) = g-^^{p\^l,^) - jR{p,tj) . (9) 

where the loop function J is given by 

J(p'. A) = -■''/ tI^C-'ts Sf {/ + IH') 7B CSf (-i + in/) 

= i f' Alj^ ,10, 

•^ Jo p-^ -p'^ -le 



in terms of the relativistic nucleoli propagator Sf{p) and W^ = 4 (m^+p^) = s. 
The renormahzed loop Jii{p,fi) = m(fi + zp)/(2 7r) follows after the proper 
non-relativistic expansion of the spectral density: 



mp I p^ mp f p^ //t f P^ 



ZTT V m^ 2tt \ 2m^ \m^ J J 

It is important to observe that for p < A the renormahzed coupling function 
can be expanded in powers of the momentum p with 

m f^ dp^ p m [°° dp^p p^ + ^x^ 



g~^\p\K-^,)^g-\p^)-^^j^ 



p^ + /i^ 27r2 J f^2 p^ — p^ p^ + /i^ 





2 / / J2 



= »-V)-^(A-i.)+^^(l + C.^i,jj (12) 

where we dropped 1/ra correction terms. Our result (|l2| ) demonstrates that all 
cutoff dependence can safely be absorbed in the bare coupling g(p) provided 
that it is legitimate to identify the typical cutoff A with the natural scale of 
the bare coupling function. 



Equation (12) exhibits an interesting effect. Starting with a perfectly 
natural coupling function gljP') it potentially generates an anomalously large 
renormahzed coupling function gi^{p^,A;^), however, only if the subtraction 
scale /x ~ m^ is taken to be small. Then for an attractive natural coupling 
g < the g~^ term is potentially canceled by the to A term such that the 
renormahzed coupling gji will be anomalously large. We stress that this is 
a familiar phenomenon since it is long understood that an attractive system 
may dynamically generate 'new' scales which are anomalously small (like for 
example the binding energy of a shallow bound state). What is intriguing here 



is the transparent ancLsipple mechanism provided by (12) to generate such a 
'small' scale naturally BqO. 

We close this section with a few simple examples demonstrating how to 
expand diagrams induced by the relativistic chiral Lagrangian and correct an 
erroneous statement made during my talk. It is sufficient to discuss scalar dia- 
grams since any diagram with Dirac structure can always be reduced to scalar 
loop function by simple algebra. Consider for example the scalar master loop 
functions which one naturally encounters in the evaluation of the relativistic 
pion nucleon box diagram: 

h,N{s) = I J -j^ S^{1 + i Q) Sn{P - I) Sn{K + I) , 
/4,jv(s, t) = -i f -^ S^{1 + \Q) SS -\Q) Sn{P - I) Sn{K + /)(13) 



with Si{p) = 1/(P^ ^ "^f + ie)- For technical convenience we introduce the 
four vectors P, K and Q as fohows: 

with (P + X)2 = s and Q^ = t hut P ■ Q ^ K ■ Q = 0. The master functions 
are easily calculated by means of dispersion techniques: 

J4m2 71" S' ~ S 
h^N[s,t) ^ / ; . (15) 



4 m^ 



TT S' — S 



The finite cutoff A < m is kept at this intermediate stage (it will be removed 
later) in order to mathematically justify the l/m expansion. We derive the 
s-channel spectral densities: 

1 1 1 / m'i — u — t 

P3,n{s) = — . „ , In 

Ibn y/4m^ - u-t \/-u - t 



, , 11 1 1 , t(u + t) , , 

P4jv(s,i) = - 1= , I — artanh < / -— — (16) 

'^^ 4^ V^ V4m2 - u - t V^Ml V ^(^'^) 

with M + t = 4 m? — s and 

6(s,t) =4TO^-4m^(u + t) + f (u + t) . (17) 

The s-channel spectral densities are expanded in powers of the inverse nucleon 
mass. At leading order we then derive: 

1 1 / (:/ZlEI\ 1] (\ " + ^ 

167rm ^J—u — t I \ m^r / 2 \ m^ 



. / N 1 1 / 1 / tiu + t) 

h N(s,t) — : I artanh W -— — 

87rm ^-i6(s,i)l I V b{s,t) 



1 

— arctan 
2 



1 

— arctan 
2 



^-tb{s,t) _ y/-U-t / ^ 

4 ?Ti^ TTItt V 4 TO 

4toJ TOtt V Am^J I I 



where we expand in K/m^ u/m and t/m. Note that the (I/ttt,)" correction 
terms of Iz^n{s) and I^^Nis, t) can be evaluated in closed form (up to a scheme 



dependent subtraction polynomial) by simply replacing the 1/m factor in (18) 
by the factor 2/y/s. This is an immediate consequence of the representation 
( |l5| ) . We find this example highly instructive since it clearly demonstrates how 
the difference of the full relativistic loop function and the l/m expanded loop 
is lumped efficiently into the subtraction polynomial. 

It is instructive to confront our relativistic scheme with the potential ap- 
proach. Inherent of the potential approach is the use of the static pion ex- 
change. The pion propagator is expanded 

SAloJ) ^ 7j-^ - ,„ J° ,,, +0 (/^) (19) 

where it is commonly argued that terms involving the energy transfer l^ are 
suppressed by l/m. Mathematically it is not quite obvious that such an ex- 
pansion is justified and in accordance with covariance since the expansion is 
highly divergent due to multiple powers of the energy transfer ^o- We evaluate 
the master pion nucleon box function /4^Ar(u,t) in the limit of static pions 

r(static)f^ -^rs I f (Pi 1 1 1 

Ian {P.P)^TZ 



(20) 

where we chose the center of mass frame with p^ = p'^ and p-ff — p^ cos and 
s = 4 (m^ +P^)- Note that we dropped terms suppressed by 1/m^. Explicit 
evaluation of the integral (120) leads in fact to the same expression derived 
before in (p^). It remains to be shown that also the (l/m)" correction terms 
evaluated in such a scheme agree with those of the relativistic approach. Note 
that in any case it becomes more and more tedious to evaluate the (I/to)" 
correction terms using (p^). 

In the algebraic reduction of the pion nucleon box diagram one encounters 
a further master loop function 

hAt) = * / 7^ s^ii + IQ) sAi -^Q) Sn{i + K) (21) 

also relevant for the pion-nucleon scattering process. The triangle loop, /3.7r(t), 
is finite and can be evaluated by means of a dispersion integral: 

.3,.(,). /"""■""' i^^ „2, 



with the t-channel spectral density: 



P3,.(i) = — ,^,^ ^ ;, arctan ^Am^^ \ / (23) 

for i > Arn^. Naively one may consider the spectral density in the heavy 
nucleon mass limit with 

P3,.W = — ^e(t-4m2) (24) 

and derive 

'^-^^)=167;;^V^"'^^^""(2^) ^25) 

for t < 0. We point out, however, that the 1/m expansion of the loop integral 
l3,Tr{t) is subtle. This can be anticipated from the fact that the 1/m expansion 
of the t-channel spectral density leads, at subleading orders, to singularities in 
the dispersion integral at i' = 4 mjR. Note that the 'naive' leading order result 
( Pq ) is in agreemeiiLjidth the corresponding expression derived in the heavy 
baryon framework 1^^"^ An alternative way to perform the non-relativistic 
expansion is through the exact identity 

h,^ (t) = /3, JV (4 m^ - i) - (m2 - i t) h^N i4m^-t,t) (26) 

which holds for the full relativistic loop functions. It expresses the desired loop 
integral Is^nit), in terms of the master loop functions /3.jv(s) and Ii^N{s,t) 
which posses a well defined 1/m expansion. The non-relativistic expansion of 
l3,Nis) and /4,Ar(s,t) then confirms (|25|). We suggest that the proper (l/m)" 
correction terms of h^Trit) arc in fact accessed via (p6|). In terms of the t- 
channel dispersion relation (23) this amounts to only including the l/m cor- 
rection terms from the factor l/\^Am? — t but use the large nucleon mass limit 
of the arctan-term in (p3|). 

The few examples presented here clearly demonstrate the usefulness of the 
relativistic chiral Lagrangian. The required 1/m expansion of Feynman di- 
agrams is straightforward and in fact performed most economically applying 
dispersion relation techniques introduced by Cutkoskytj. These techniques 
keep covariance manifestly andpaxe applied easily also to higher loop diagrams 
involving pion production cuti2l. Note that dealing with relativistic Feyn- 
man diagrams avoids the cumbersome splitting of Jiions into static, soft and 
radiation pions as suggested by Mehen and Stewarcj. 

'^Note that such singularities disappear in the chiral limit with m-n = 0. 

10 



3 Chiral expansion scheme for nuclear matter 

An attempt to apply the generalized chiral power counting rules oiuu to the 
nuclear many body problem quickly reveals that even though the pion dynam- 
ics remains perturbative the local two-nucleon interaction requires extensive 
summations. This is an immediate consequence of the chiral power given to 
the renormalized coupling Qr ^ Q~^ ■ In nuclear matter at small density with 
kp — 100 MeV and p = 2 k'p / {?> ■k'^) it is not sufficient to sum the particle- 
particle ladder diagrams of the local two-nucleon interaction. We find that 
one must sum simultaneously also the particle-particle and hole-hole ladder 
diagrams including all interference terms where a self consistently dressed nu- 
cleon propagator must be used. According to the generalized counting rules 
the pions can then be evaluated perturbatively on top of the above described 
parquet resummation for the local two-nucleon interaction^. In terms of vac- 
uum scales this leads to an expansion of the energy per particle, Eikp), of the 
form 

with typical small scales, Ag, like y/m^ ejj where cd ~ 2 MeV is the deuteron 
binding energy, and typical large scales Kl — ^tt f-K — 1 GeV. The expan- 
sion coefficients En are complicated and hitherto unknown functions of the 
Fermi momentum kp. They can be computed in terms of the fcee space chiral 
Lagrangian furnished with a systematic summation techniqueE^. 

In this work we pursue a somewhat less microscopic approach in spirit 
close to the Brucckner scheme but more systematic in the sense of effective 
field theory. Since the typical small scale As is much smaller than the Fermi 
momentum kpfi — 265 MeV at nuclear saturation density, one may expand 
the coefficient functions £"„ around kp.o in the following manner 



£4^.i^^B.fe¥^vgfi!^'(ie) (ij-^i ^ (-) 



m^ ^s J \m^ A5 



Note that we do not expand in the ratio mT^jkp. If one expandjEd also in 
this ratio raT^lkp one would arrive at the Skyrme phenomenology E3 applied 



Our expectation that the nuclear equation of state can be calculated microscopically solv- 
ing the parquet theory for the local interaction and then include pions perturbatively may 
be somewhat too naive. A refined formulation which incorporates iselevant subthreshold 
singularities of the vacuum nucleon-nucleon scattering amplitude (sccq) may lead to a more 
involved treatment of pionic effects. 

11 





Figure 1: Leading contribution of chiral order Q°. 

successfully to nuclear physics many years ago. It should be clear that this 
scheme is constructed to work around nuclear saturation density but will fail 
at small density. We note that also conventional approaches like the Walecka 
mean field or the Brueckner scheme are known to be incorrect at small density. 
In terms of diagrams one may arrive at the proposed expansion (|2g) by a 
perturbative evaluation of the pionic nuclear many-body effects on top of the 
the high density limit of the parquet resummation for the local two-nucleon 
interaction. The underlying assumption is of course that the typical small scale 
in the parquet, Ag, is such that the high density limit is reached already at 
densities much smaller than the nuclear saturation density kp^Q. Technically 
our scheme can be generated by the effective Lagrangian density 

C^ntikp) = If iV75 r ■ {d^. 7?) • tN 

+ ^ (5o(fcF) + \ 5i) (iV75 T2 C-1 N') (iV* Ct2 75 N) 
+ ^ (5i(fcF) + Jsi) (iV7^fr2C-i7V*) (iV*CT2f 7^ iV)(29) 
where the couplings go — goikp) and gi — giikp) are density dependent. A 



ikp) 



M) applying suitable 



more systematic derivation of the expansion (^7|) at 
resummation techniques will be presented elsewhereE2 . 

The chiral power counting rules are simplified significantly as compared to 
a fully microscopic scheme. The presence of a further small scale kp ^ Q ^ m-„ 
does not anymore generate an infinite tower of diagrams, to be considered at 
a given chiral order, since by construction the troublesome local 2-nucleon 
vertex need not be iterated, i.e. terms proportional to g^ikp) and gi{kp) 
with n > 1 are already included in goikp) and gi{kp) and therefore must not 
be considered. Note that here we count go,i ^ Q° since the non-perturbative 
structures like the deuteron, which give rise to the anomalous power go,i ~ 
Q~^ in the vacuum, are long dissolved at densities close to nuclear saturation. 
The pion dynamics, if properly renormalized, remains perturbative like in the 
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Figure 2: Subleading contribution of chiral order Q 



vacuum case. Consider for example the two foop diagrams depicted in Fig. || 
where the nucleon Hnc with a 'cross' represents the projector onto the Fermi 
sphere 

A^^^(p)-(7•p + m^,)2 7r^e(po)<5(p'-m^)e(4-p2) . (30) 

The first diagram in Fig. |l| is proportional to go,i{kF) kp and is therefore 
ascribed the chiral order Q^ since the effective vertex go^iikp) ^ Q^ carries 
chiral power zero. The second diagram, the one pion exchange contribution, 
is also of chiral order Q^ since it is proportional to k^p multiplied with some 
dimensionless function f{m-,^/kp). In Fig. g we collected all diagrams of chiral 
order Q^ . Here wc introduced two types of 2-nucleon vertices. The filled circle 
represents the full vertex of (29) proportional to 5o,i +5^/4 and the open circle 
the counter term proportional to 5^/4. The dashed line is the pion propagator 
and the directed solid line the free space nucleon propagator. We point out 
that the diagrams b), d), f), h) and j) in Fig. || are divergent. The leading 
chiral contribution of the sum of all diagrams, however, is finite after including 
the appropriate counter term in the spin triplet channel. 

It is instructive also to display the non relativistic form of the Lagrangian 
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Figure 3: The equation of state for isospin symmetric nuclear matter. 

density (129) as suggested by Weinberg EJ in the context of the nuclcon-micleon 
scattering problem 



A„t(fci.) = |^iV 



N 



(a • V Vtt • t] 

+ i (90{kF) + ^gyj {N ff (T2 T2 N'){N't2 (72 a N) 

+ Yp (■9l(^^) + i5A) (^^^2 TT2N') {N'T2f<J2N) (31) 

with N now a two component spinor field and N — N'^ for notational conve- 
nience. According to Weinberg the two-particle reducible diagrams are to be 
evaluated with static pions. 

We therefore evaluate all diagrams of Fig. || with the appropriate non 
relativistic interaction vertices and static pion propagators. The solid line 
with a 'cross' now represents the non relativistic limit of (^0|) and the directed 
solid lin£_the free non relativistic nucleon propagator. For technical details we 
refer too. Both schemes indeed lead to identical results for all diagramsy. 

In Fig. H our result for the isospin symmetric nuclear equation of state is 
shown. The relevant coupling go + gi is adjusted to obtain nuclear saturation 

0. 



'^Note that here we correct an error in 



14 



at kpfi = 265 MeV. We emphasize that the coupUng functions go.iikp) are to 
be determined from the nuclear equation of state. However there is a strong 
consistency constraint: according to our scale argument ( P8|) the density de- 
pendence of the couplings go,i{kF) rnust be weak for kp larger than typical 
small scales integrated out. If the nuclear saturation required a strong density 
dependence our scheme had to be rejected. The density independent set of 
parameters go + gi — 3.2, gA — 1.26, TOtt — 140 MeV and /^r — 93 MeV give 
an excellent result for the equation of state. The empirical saturation density 
kpfi — 265 MeV and the empirical binding energy of 16 MeV are reproduced. 
The incompressibility with, k ~ 215 MeV is also compatible with the empirical 
value (210 ± 30) MeV ofEJ. 

4 The chiral order parameter (qq) 

The quark condensate, {qq){p), is an object of utmost interest. It measures 
the degree of chiral symmetry restoratkm in nuclear matter. Furthermore it is 
lEB important input for QCD sum rulesLj or the Brown- Rho scaling hypothesis 
c3. According to the Feynman-Hellman theorem the quark condensate can be 
extracted unambiguously from the total energy, E{p), of nuclear matter once 
the current quark mass dependence of E{p, mq) is known 

{qq){p)^^.j^E{p,mQ) (32) 

V oniQ 

with niQ = TTiu — md- Recall that E{p)/V — {mN + E{p)) p is determined 
by the nuclear equation of state E{p). Since our chiral approach was set up 
to treat the pion dynamics and therewith the current quark mass dependence 
of the equation of state properly it is very much tailored to be applied to the 
quark condensate. 

It is convenient to consider the relative change of the quark condensate 
since it is renormalization group invariant: 

(gg)(p) ^ -, _ ^ttN P _ a^{p)p .„„. 

{qqm ml,n 2 m. f^ ^ ^ 



The second term in ( p3[ ) follows from the nuclcon rest mass contribution to the 
total energy of nuclear matter together with the definition of the pion nucleon 
sigma term, S^at: 

S.JV = TTiQ {N\qq\N) - niQ —^ . (34) 
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Figure 4: The quark condensate in isospin symmetric nuclear matter. 



The last term in (^3|) measures the sensitivity of the nuclear equation of state, 
E{p,mTr), on the pion mass 



a^{p) 



{qq){0) dmq 



E{p,m^) = {l + 0{ml))- E{p,m^) (35) 



drriTr 



where we consider now the current quark mass mg = mQ^rriT^) as a function 
of the physical pprpiiass. We emphasize that the second term in ( p3|) writ- 
ten down first inEjOo does not probe the nuclear many body problem and 
therefore should be considered with great caution. It is far from obvious that 
this term is the most important one at nuclear saturation density. The only 
term susceptible to nuclear dynamics is the last term in (33). 

In Fig. y we present our result for the quark condensate in nuclear matter. 
We confront the 'leading' term driven by 'S-nN — 45 MeV with our result. 
The inclusion of pionic many body effects leads to a significantly less reduced 
quark condensate in nuclear matter. Note that here we do not include an 
explicit r?T,^ dependence of the coupling go^i{kp) simply because there is not 
yet any reliable estimate available. Our resultaxonfirm calculations performed 
within the BruecknerLj and Dirac-BruecknerLiJ approach qualitatively insofar 
that the nuclear many body system appears to react against chiral symmetry 
restoration. They have strong implications for the QCD sum rule approach of 
hadron properties in nuclear matter. 
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